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We investigate the decay of metastable states with broken CP-symmetry which have 
\ recently been proposed by Kharzeev, Pisarski and Tytgat to form in hot hadronic 

matter. We consider the efficiency of the amplification of the r/'-field via parametric 
resonance, taking the backreaction into account. For times of the order t ~ 10/m, we 
find a particle density of about 0.7 //m 3 and a correlation length of £ ma;c ~ 2.5 fm. 
' The corresponding momentum spectra show a non-thermal behaviour. 
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1 Introduction 



Recently, Kharzeev, Pisarski and Tytgat [1,2] presented the following idea, 
based on the nontrivial topological structure of the QCD- vacuum [3]: Metasta- 
ble states which act like regions with a non-vanishing QCD vacuum angle 
6 may be excited when matter undergoes the deconfining phase transition, 
provided it is of second order. In these false vacua, parity (and also CP) is 
spontaneously broken, a quality that could lead to experimental signatures 
like parity odd correlations of the produced particles. 

In the present work, we investigate the production of ?/-particles during the 
decay of the CP-odd metastable states, concentrating on the amplification of 
the low momentum modes by parametric resonance when the zero mode rolls 
down from the false to the true vacuum. This mechanism plays an important 
role for particle production in inflationary cosmology [4-11] and has also been 
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investigated for the formation of Disoriented Chiral Condensates [12-15] and 
for Axions [16]. Here we discuss the efficiency of the amplification mechanism 
for the production of 77 from the decaying metastable bubble. As a main re- 
sult we derive the corresponding momentum spectrum showing non-thermal 
behaviour. The integrated density of particles produced is estimated to be 
0.7/m" 3 . 



2 Metastable states in hot hadronic matter 



The idea is based on the effective Lagrangian of the Witten-DiVecchia-Veneziano 
model [17-19] 
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eff 



f 2 
T 
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(1) 



which describes the low-energy dynamics of the pseudoscalar mesons in the 
large A^-limit of QCD. In the following, we consider the "real world"-vacuum 
angle = 0. Formation and decay of a non-zero ^-vacuum and its signatures 
are investigated in [20]. The Nf x iVf-matrix U in Eq. (1) describes the meson 
fields, and the term containing the mass matrix M provides the explicit soft 
breaking of chiral symmetry due to the quark masses. M can be written as 



M, 

t4 



pqSiji where the diagonal elements are chosen as fxf 



f4 



mz 



and 



2m\ 



m 2 .. With the parametrization U = exp(i(J) / f w ) , the matrix 
representing the singlet and the octet meson fields reads 



:?7il + 



V2K~ 



V2K° 



(2) 



The last term in the effective Lagrangian reflects the £/4(l)-anomaly by giving 
a mass to the singlet, even in the chiral limit of vanishing quark masses. The 
parameter a = 2Nf\ Y M/fl represents the topological susceptibility. 

Finite temperature is introduced into the model via the temperature depen- 
dence of the parameters: For T = 0, the value a = m 2 + mL — 2m 2 K ~ 
0.726 GeV 2 can be determined from experimental results [17-19], and fi 2 = 
(fif +1^1 + /4)/3 = (ml + 2m 2 K )/3 ~ 0.171 GeV 2 . The corresponding value of 
the pion decay constant is f w ~ 93 MeV. 
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When the temperature approaches the temperature Td of the phase transition, 
a(T) goes to zero, indicating the effective restoration of the L^(l)-symmetry 
for T > T d with enhanced rj - and rj' - production as a possible signature 
[21,22]. We use the temperature of the deconfinement phase transition, T d) 
following the assumption of [1] that any other phase transition of this model 
also occurs at According to mean field estimates [1], a(T) oc (Td — T) and 
/x 2 (T) oc (Td — T)~ l l 2 near the phase transition. 

In order to get an analytically tractable approximation that allows some in- 
sight into the physical processes, we consider only the singlet, which is the 
main component of the 7/-meson. Using 7V C = 3, the effective Lagrangian 
simplifies to 

£ s = l^mWm) + \flf cos (yf |] - a - v l (3) 



The singlet effective potential 



t/ I 'I \ - (^\ r] a rf 



where we have dropped the index 1 in the notation of the singlet and intro- 
duced / = ^3/2/V for convenience, can be considered as a projection of the 
full effective potential. It takes three qualitatively different shapes, depend- 
ing on the temperature, more specifically, on the value of a(T)/fi 2 (T). For 
large values of a//! 2 , corresponding to small T, the potential is a more or less 
deformed parabola (see Fig. la). For a/jj? = 0, i.e. T > T d , the potential 
is periodic with minima at rj/f = 0, 2n, An..., which are all equivalent true 
vacua (see Fig. Id). And for a small, nonzero value of a/fi 2 , corresponding 
to a temperature just below the phase transition, metastable states appear 
in the potential which are distinct from the true vacuum (see Fig. lc). One 
recognizes that the false vacua are odd under parity transformation, but even 
under charge conjugation, and therefore odd under CP [1]. 

One can show that the effective potential exhibits metastable states also in the 
case of nonvanishing 7i—,r] 8 — and K— fields, provided that m n ^ [23]. But 
only the singlet is responsible for the formation of these states, in the sense 
that it is the only field which appears in the term of the effective potential 
which includes the topological susceptibility. 



3 The decay of the metastable states 
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3. 1 Dynamics of the model 



Applying the model to the description of a heavy ion collision, we suppose that 
during the phase transition, the potential changes its shape from the parabola 
(Fig. la), corresponding to the low temperature phase, to that of the high 
temperature phase (Fig. Id), and back to the low temperature phase, but 
this time including metastable states that have been created at a temperature 
just below the phase transition. After the phase transition, but still at high 
temperature, the 77'— field may be trapped in a metastable state, forming a 
"CP-odd bubble inside the hadronic phase" in the language of [1]. Considering 
the effective potential as a function of temperature, one recognizes that there 
exists a value a/ n 2 {T sp ) for which the local minimum turns into a saddle point 
(see Fig. lb). For the lowest false vacuum, this happens for a/fx 2 (T sp ) = 0.217 
at the value (rj/f) sp = 4.493. Assuming the proportionality a(T)/fi 2 (T) oc 
(Td — T) 3//2 to hold even far from T^, the temperature corresponding to this 
special case of the potential can be estimated to be T sp ps 0.86 Tj. Once 
the saddle point temperature is reached and the barrier has disappeared, the 
rj'— field starts rolling down into the true vacuum and oscillating around it, 
while energy is transferred from the mean value of the field into its quantum 
fluctuations. The resulting growth of the occupation numbers of the quantum 
fluctuations is interpreted as particle production with respect to the true, CP 
even, vacuum, whereas the zero mode is damped. This process is investigated 
in more detail in the following. 

A crucial assumption for this picture is that at least after the saddle point 
temperature has been reached, the effective potential changes slowly compared 
to the motion of the field, so that the oscillations can be considered as taking 
place in a static potential, frozen at T = T sp . Accordingly, the ratio a//i 2 = 
0.217 is kept fixed during the evolution in time, and the value of / as well. 
The solutions we find (see below) show that the zero mode reaches the true 
vacuum for the first time after 1 — 2 fm, whereas the relevant time scale to 
be considered for the particle production is of the order of 10 fm. 

3.2 Evolution equations 

To study the dynamics of the decay process, we decompose the field into its 
expectation value ip(t) = (r)(x,t)) and fluctuations x about it, 



with (x(x,t)) = 0, so that the "physical" fluctuations are given by the usual 
expression ((% — (x)) 2 ) = (x 2 )- With this decomposition, following the stan- 
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dard method as described in, e.g., [5,6,8,11], we obtain two evolution equations 
from the effective Lagrangian (3). The effect of the quantum fluctuations at 
lowest order is taken into account with a Hartree-type approximation [5,8,11] 
that consists of the factorization x 3 ~* 3(x 2 )x an d X 2 ~^ (x 2 ) an d the self- 
consistency condition represented by Eq. (7) given below; terms of higher order 
in x are neglected. Taking the expectation value (), we arrive at the equation 
for the zero mode 



written for the dimensionless mean field f/f, where the dimensionless time 
variable r = fit is introduced. The expectation value of the fluctuations 



is calculated from the mode functions Xfc( r ) which are known from the Fourier 
expansion of the field. The integration is done over the interval of momenta 
for which the quantum fluctuations are enhanced. 

We use the definition (7) for the back reaction, although it does not fulfill 
(x 2 )o — 0; because it reflects the background of the calculation, the closed- 
time-path technique which is usually applied for the description of out-of- 
equilibrium dynamics in field theory. In this framework, the fluctuations are 
defined as the two-point correlation functions or equal-time Green's functions 







[5] 





The mode functions satisfy the mode equation 



dr 2 



Xk(t)+u 2 k (t)xk(t)=0, 



(9) 



written for the dimensionless mode functions Xk = y^Xfc, where 




(10) 
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is the time-dependent frequency squared with k = k/fi. They are used for the 
calculation of the particle production: The induced particle number density 



n(r) = / 



d 3 k 



n k (r) 



(11) 



(2tt) 3 



is obtained from the spectral particle number density 



n K (r) 





)l 



2 



) 



2 



1 



(12) 



which can be calculated from the mode functions. Since the unstable solutions 
of the mode equations contribute to the resonant particle production, the 
integration in (11) is carried out over the momenta corresponding to these 
amplified solutions, as in Eq. (7). For practical purposes, we use ur K = /t 2 +a//x 2 
in definition (12), because it keeps the particle number positive definite for all 
values of r and minimizes the fluctuations in the particle number. Accordingly, 
the initial values of the mode functions are chosen as x«(0) = 1 / and 



Xk(0) — — iy& K /2 such that n K (0) = 0. However, the results depend only 
weakly on the choice of oo K both in the initial conditions and in the definition 
of the particle number. The choice for the initial conditions of the zero mode is 
suggested by the model: is "released" right below the saddle point (p sp / f = 
4.493 with zero velocity, i.e. ip(0) = 0. For the numerical results shown below, 
we used the starting value <^(0)/ / = 4.2. 

3.3 Solutions and results 

Considering the coupled system (6) and (9) with (10), one recognizes that the 
expectation value transfers energy to the mode functions via a time dependent 
frequency, and the mode functions in turn modify the equation of motion for 
the zero mode. A natural first step in the solution of this system is to neglect 
the fluctuation terms in both equations, obtaining the "classical" evolution 
equations with (x 2 )r = where the mode functions do not react back on the 
zero mode. If the variable frequency 



in the mode equation depends periodically on time, according to Floquet's 
theorem [24] it has quasiperiodic solutions of the form 





(13) 




(14) 
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where the characteristic exponent \i K depends on u Kt d, and P(t) is a periodic 
function of the same period as ou KtC i, usually normalized to have unit amplitude. 
For certain values of the parameter oj k>c i, there exist so called resonance bands 
that contain the exponentially growing (or resonant) solutions with a real 
characteristic exponent /i K . These resonance bands are known, for example, 
for the solutions of the Mathieu- or the Lame-equation. 

An analytic solution of the classical system can be found with a Sine-Gordon- 
equation [7], which consists essentially of setting a//i 2 in our case. The 
corresponding approximate zero mode equation is given by 



where p = 1.43 represents the radius of the potential and is needed to adjust 
it to the original potential. It has the solution 



with e = l — cos(tp(0)/pf), and m = e/2 is the square of the modulus of the 
Jacobian elliptic function cn. The resulting mode equation is a Lame-equation 
with the resonance band [7] < k 2 < e/2. The resulting fj, K is shown in Fig. 2 
in comparison with the result of the numerical calculation of \i K for the solution 
of the classical mode equation. In the range of the single resonance band of 
the Lame-equation, < k 2 < 0.97, the numerical calculation yields three 
resonance bands, and at least two for k 2 > 1. (This structure reminds of the 
Mathieu-equation which has more than one resonance band.) The maximum 
value (J, K;m ,ax — 0.42 is found around k 2 ~ 0.1, and the maximum value for the 
Sine-Gordon-solution, p s K 9 max — 0.5, yields an acceptable approximation . 

The momentum spectrum of the produced particles reflects the structure of 
the resonance bands. When the general expression for the mode functions (14) 
is inserted into Eq. (12), it is obvious that the particle production is essentially 
determined by /i K : 



In the classical approximation discussed above, the back reaction of the fluctu- 
ations on the zero mode (and also on the fluctuations themselves) is ignored. 
This implies that the approximation is not energy conserving, which can im- 
mediately be seen from the zero mode solution which is not damped although 
it is supposed to transfer its energy to the mode functions, which indeed grow 




(15) 




(16) 




(17) 
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exponentially. Since the mode functions are not decoupled from the zero mode, 
at least the other direction of energy transfer is reflected by the approxima- 
tion: The larger the initial value tp(0), the more energy can be transferred, 
and the larger the mode functions and the particle number grow. 

The evolution equations including the back reaction are solved numerically. 
For the evaluation, we use f(T sp ) ~ /(0) and n(T sp ) ~ 676 MeV. Because 
the value of // in our approximation is only fixed by the value of T sp , which 
depends on our restriction to the singlet case and may not correspond to the 
real physical situation, we also try /x(T = 0) = 412 MeV as the other limit of 
the range of possible values for fi(T). 

Although the inclusion of the back reaction destroys the parametric resonance, 
it does not suppress the particle production completely. The particle number 
(Fig. 3) reaches its asymptotic value for t en d ~ 10 fm. At about the same time, 
the value of the back reaction term (x 2 )r/2/ 2 grows larger than 1, indicating 
that the approximations made become inappropriate for the description of the 
dynamics for larger times (cf. Eq. (6)). One recognizes from Fig. 3 that the 
choice of /x(T) has some influence on t end ; the asymptotic number density of the 
produced particles lies between 0.7 fm~ 3 and 0.9/m" 3 . The absolute number of 
particles can be calculated from this density if the size of the CP-odd bubbles 
is known. For example, for a domain with radius 5 /to, as suggested in [25], 
we estimate a yield of 90-100 particles from our model. 

A measure for the size of correlated domains of particles is the correlation 
length £, which can be calculated from the two-point correlation function 



from which we estimate the correlation length for t en d- As a lower limit for the 
domain radius, we obtain £ ~ 2.4 fm for /x(T = 0) and £ ~ 1.3 fm for fi(T sp ). 

The momentum spectrum of the produced particles at t en d ~ 10 fm is shown 
in Fig. 4. Its maximum is located at about the value of k as the maximum of /j, k 
in Fig. 2, which had been calculated for the system without back reaction. The 
spectrum deviates significantly from a thermal Bose-Einstein distribution. 




(18) 



where x — jir denotes the dimensionless length. For large x, 




(19) 
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Figure Captions 



Fig. 1. Form of the singlet effective potential for different temperatures: (a) T = 0, 
(b) T = T sp , (c) T sp <T< T d , (d) T d <T . 

Fig. 2. The characteristic exponent \i K for the solutions of the classi- 
cal evolution equations (solid line) and for the analytical solution of the 
Sine-Gordon-approximation (dotted line). 

Fig. 3. Particle number density per /m 3 as a function of time for 
fi(T = 0) = 412 MeV (dashed line) and fi(T sp ) ~ 676 MeV (full line). 

Fig. 4. Momentum spectrum of produced particles at t end ~ 10 fm, including back 
reaction and using /j,(T = 0) = 412 MeV. 
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